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I. INTRODUCTION 



Typically, given a model of an interacting spin system 
on a (i-dimensional regular lattice, the first approach in 
understanding its properties is through the Weiss mean- 
field theory. Equivalently, one can neglect the underlying 
lattice structure completely, where each spin interacts 
with a finite number of particles, and assume that all 
spins interact with each other. In other words, one sub- 
stitutes the original model by its (exactly solvable) fully- 
connected version. This crude step may fail either com- 
pletely, for instance predicting phase transitions where 
there are none, or give a rather good qualitative picture 
of the physics of the model. In order to improve results 
both qualitatively and quantitatively one must then use 
sophisticated techniques, such as renormalization group. 
During the last two decades, since the seminal paper 
of Viana and Bray Q, we have witnessed that it is 
still possible to consider interesting exactly solvable ver- 
sions of the original model, other than the fully con- 
nected ones, while keeping the coordination number fi- 
nite. Since then, interacting spin systems on random 
raphs have attracted quite some attention in spin glasses 
Jill El El IE pi neural networks 0)11,111 an d small- world 
networks fialll|. 
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Surprisingly, during the last decade we have encountered 
interesting problems that can be mathematically formal- 
ized precisely as interacting spin systems on random 
graphs. Examples can be found in optimization prob- 
lems H2, [H, [3, lilSll. error-correcting codes [H, H3, 
and cryptography |l9i l2(i| . 

In order to tackle this type of problems the replica ap- 
proach and cavity method, fairly well understood for fully 
connected systems pi], had to be reconsidered 0, E|: 
while in the one pure-state approximation plj a few or- 
der parameters suffice to describe the thermodynamics in 
the fully connected case, this no longer holds as soon as 
the coordination number is finite. In the latter case an 
infinite set of p arameters, or a function order parameter, 
is required p2 |. 

Unfortunately, while most analytic studies of interact- 
ing spin systems have considered spins of a discrete na- 
ture, e.g Ising spins, much less attention has been given 
to systems with real-valued variables and in particular 
their model realization on random graphs. Needless to 
say, such systems are not only more realistic from a fun- 
damental physical point of view, but are also of special 
relevance in several contexts, as for instance in Josephson 
junctions arrays piil l24l |25| . Kuramoto models pa. |27|. 
etc. 

The goal of this paper is then twofold. Firstly, we wish 
to extend the cavity method to spin systems with real- 
valued variables. Secondly, and following the steps in 
p8| . we would like to study here the thermodynamical 
properties of Kuramoto models in small-world networks 
|29|. Small worlds, more than just curious structures, are 
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apparently universal and can be observed in many differ- 
ent circumstances: from linguistic and social networks 
to the world- wide- web (for a recent review see e.g. |3Cf ). 
A large body of work has been devoted to the study of 
small-world networks, mainly numerically (see e.g. |3l|) 
and to a lesser extent analytically [lfl Hl| . 
This paper is organized as follows: in Section [H] we in- 
troduce the XY spin model in small world networks. To 
solve it we generalize the cavity method in Section IIIII 
to real-valued dynamical variables to specific graph in- 
stances and in the ensemble case. We then use the cavity 
method on the XY small-world spin model in Sect ion ITVl 
we first apply bifurcation analysis to calculate the con- 
tinuous transition lines exactly. We also use population 
dynamics and compare the results with simulations. We 
finish Section IIVI by studying the proliferation of pure 
states. The last section is for conclusions. 



II. DEFINITIONS 



We study a small-world network of N interacting XY 
spins <Tj = (cr?, of) with i = 1, . . . , N and with ||<Ti|| 2 = 
<rf = 1. At equilibrium the system interacts according 
to the Hamiltonian 



H({*i}) = -J J2<n 



w E < 

(i,j')es 



•U 



(1) 



The extension of the cavity method to any type of dy- 
namical variable was worked out by one of us in |32| and 
we would like to recall it here. To do so, let us first 
consider a general version of the Hamiltonian , viz. 



H({<Ti}) 



J E ■ 



U,: 



(2) 



where the nature of the spins is unspecified (Ising, real 
variables, spherical, multicomponent, and so on) and 
we will require only that the partition function is well- 
defined. Let us assume that the graph Q is tree-like (short 
loops are rare). It is then clear that the spins belonging 
to the neighborhood of i are mainly correlated through 
the spin at site i. If, however, we consider the fictitious 
situation of removing the spin at site i, then in the re- 
sulting cavity graph the joint probability distribution of 
its neighbors approximately factorizes and we can write 



where the first term represents a ferromagnetic interac- 
tion along a one dimensional chain while the second one 
describes a chiral interaction. The factors are or- 
thogonal matrices drawn from some general distribution 
-D(U). Here Q = (V,£) represents a graph with vertex 
set V = {l,...,N} and edges £ 9 £ C (V, V). We 

denote by di the neighbors of site i and for any subset of 
vertices A C V we write a a = {o~i\i S ^4}- 



III. GENERALIZED CAVITY METHOD 
A. Cavity method on an instance 

The cavity approach for diluted disordered systems 
with Ising-type variables has been developed in @, 0- 



Bct.hc 
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(3) 



where Pj (iTj) denotes the cavity probability distribution 
of finding spin j in state <jj in the absence of site i. This 
is usually called the Bethe approximation. Let us now 
consider the reverse situation where we start from a cav- 
ity graph without the spin at site i and reincorporate the 
spin i and reconnect all sites di to i. Then the physical 
(marginal) probability distribution for cr; on the original 
graph Q is given by 



Pi{<Ti) = y 



da-Qi e 



, (Ti-Uii rr 



; (Ti-Uie (J, 



ted- 



(i) 



[<Tt 



(4) 



r 



with Zi = J d<rPi(cr). The marginal Pj((Tj) is used to 
calculate the physical one-site quantities in the original 
graph Q. Unfortunately the cavity distributions P £ (cri) 



are still not known. To overcome this problem we notice 
that that we get closed relations for them if we reconnect 
all sites but one j £ di 
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PP{*i) = 4fi [*re«jeP J Z««*W< ^"W- Z^=fd<rP«\o) 



(5) 



iedi\j 



We see now that eqs. (JSJ) provide recursion relations 
for the cavity distributions P^\(Ti) Vi = l,...,N and 
Vj G <9i. Once these are known we can calculate any 
physical one-site quantity by using eqs. (@J. For instance 
the magnetization vector M and the spin-glass parameter 
matrix Q are given by 



i N r 

M a = jjY.J ^^(^K 



1=1 

N 



Qab = T7 X] / dcrrfT Pi{cr)P l (T)a a 
i=l ^ 



(0) 



where a, 6 labels the spin components. Note that the 
cavity distributions appearing in eqs. ||2J and JSJ) are not 
expressed in terms of parameters. We will call this choice 
of representation of the cavity equations the parameter- 
independent representation. 

In contrast, let us go now to a parameter- dependent 
representation of the cavity distributions. We denote 
a parametrization explicitly as P^{cr) — ► P(<r|/nj )■ 
Since any probability distribution is defined to be positive 
the following exponential parametrization is convenient 



(7) 



\n=0 



where fj,^ — {/4n}n>0 are ^ ne so-called cavity fields and 
{ip n (cr)} is a set of orthogonal polynomials with respect 
to a weight w(<x) 



J da u>(<r)(p n {cr)(p m {cr) 



(8) 



We will assume that the first polynomial <po((r) is a con- 
stant (the normalization constant of the probability). 
Now, using (7J in (JSJ) we find the following recursion re- 
lations for the cavity fields 

$1= E «m(M?V.U«), m>0, 
u m (/Lx|J,A) = J daw(o-)ip m (cr) (9) 

d Te /3^-AT+Er>o^ y n(T) ) 



X log 



where u(fj,f'\J, Uu) = {u m (p^\J,Uu)} rn >i are the 
messages or propagated fields. Notice that if the dynam- 
ical variables are Ising-type, then only one cavity field 
suffices to parametrize the cavity distributions and we 
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recover the method of 0, 0] . 

Let us also use the same parametrization of the marginal 
distribution Pi{<Ti) — > P(a , i\fi i ) ) where the \i i are the 
physical fields. Then from eqs. we get the following 
relation between the physical fields and the cavity fields 



£€di 



04V,U«), m>0. 



(10) 



Equations © and (|10|l are usually referred to as mes- 
sage passing equations and have the following interesting 
physical meaning: when site j S di is removed then each 
site i £ <9i\j sends a message u(/ig | J, ILj) to site i ac- 
cording to the interaction between spins I and i and the 
state of spin erg. Then the spin <7i finds its optimal con- 
figuration (through its cavity fields fif ) in the cavity 
graph by adding all messages up, eq. 0. 



B. Cavity method on the ensemble 



Looking at the expressions of the order param- 
eters we notice combinations like (l/N) X)i=i P%{cr) 
(1/N) J2iLi Pi^Pi 1 ") and so forth. It is therefore in- 
structive to introduce a density of marginal distributions 
W p hys[{f }] as follows 



W, 



phys[ 



1 N 



with §(f)("') a functional Dirac delta. Thus, we can 
rewrite the order parameters as follows 

M a = J \dP}W phys [{P}] Jd<rP{a)a a 

Q ab = J {dP}W phys [{P}] J dadr P(a)P(T)a a n . 

(12) 

Analogously, we can also define a density of cavity distri- 
butions W[{P}} 



N 



W[{P(] 



i—1 1 1 



(13) 



jedi 



If we then assume that the graph ensemble is character- 
ized by a degree distribution pk , from equations (JSJ and 
(|13fl we obtain the following self-consistency equations 
for the density W[{P}] and by using equations l@J, ltTT|) 
and (O we relate the densities W p hys[{-P}] and W[{P}] 
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fc=0 x 1 J 

oo „ 

W phys [{P}]=Y,Pk / 

7, — n •* 



•fc-i 



JJ{dP^}W[{P^]dU^(U^) 

fe 

j]{dpaw[{A}]rfu^(u,) 



'(F) 



P(<r) - - 



k-1 



11 dr e P e (r e ) 
Le=i 

k 

j I dr e Pe(Ti) 



<■:=! 



(14) 
.(15) 



Notice that the probability of choosing a cavity site ran- 
domly is proportional to the number of its bonds, i.e. 
Pkk, as reflected in the preceding equation. Finally we 
can also write self-consistency equations for the cavity 
fields fi in the ensemble. We can do this in two ways. 
Let us start from eq. Q and analogously to the defini- 
tion of the density of cavity distributions (|13fl we define 



a density of cavity fields 



N 



" i 1 1 jem 



0) 



(16) 



we then get 



J 



■fc-i 



Y\_ dfj, e w(n e )d\Je.D(lJe 
I 



U=l 



k-1 



5 [i - y^u OT (jUi|J,Ui) 



(17) 



Alternatively, we would like to get the self-consistency 
equations for w(fi) directly from those of W[{P}] as done 
in [28|. It is instructive to notice that we can get the 
transformations w(fi) <-> W[{P}] from the definition H6|) 
and of the density of cavity fields w(fj,) and cavity 
distributions W[{P}], respectively. Indeed, let us simply 



J 



assume that for a given set of cavity fields fj,[^ there is a 

unique cavity distribution P^ (a) and vice versa so that 

we are entitled to write fi^ = fj,({P^}). Then from 
eqs. lfTo)l and f|T~Sf> we obtain 



i N I 



1 N 1 r 

4 =1 1 1 3 - e 9i 



{dP}W[{P}]<5[/x-/x({P})] 

d[iw{[i)S( F) [P(cr) - P(«r|/i)] 



(18) 



respectively, which is precisely the transformation 
w(fJ,) <-> W[{P}] for switching between parameter- 
dependent and parameter-independent representations 
used in |28| . 

Finally the order parameters can be expressed as an in- 
tegral over the density of cavity fields by using the trans- 
formations l|18fl 

M a = J dfi Wp h y S (fJ,) JdaP(a\fi)a a 

Qab = J dnw phys (n) J dcrdr P(cr\n)P(T\n)(j a T b , 



where u> p hy S (/-*) is the density of physical fields. 



IV. APPLICATION TO THE XY SPIN SYSTEM 
IN SMALL WORLDS 



Let us now consider the XY model in small worlds 
defined by the Hamiltonian J2J ■ The application of the 
cavity method is straightforward. We must take into ac- 
count that the spins are XY, so that the integral over spin 
degrees of freedom is over the unit circle, J da — > f„ da 



(19) with S = {(<7 x ,a y ) e 



1}, and that the en- 
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ergetic contribution coming from the chain and from the 
graph is different. 

Indeed, let us define C(i) = {i — l,i + 1}. We will distin- 
guish whenever the neighbor j G di belongs to the imme- 
diate neighbors of the chain, j S C(i), or not, j £ di\C(i). 
Then messages transmitted along the chain between two 
spins <y and r have an energetic contribution of — Jq<t ■ t 
while for those messages along the graph the contribu- 
tion is —Jcr ■ Ut. Thus the message passing equations 
in this case simply read 

ViX = E u™(/4 l) l J o, M 
iec{i) 

+ u m {nf\J,\J it ), jedi\C{i) 

eedi\j 

$1,= E u m (fi { e i) \j ,i u ) 

eec(i)\j 

+ ^u m (/4 i) |J,U a ), jeC(i), 

iedi 



(20) 



In order to have the equations in the ensemble we define 
two densities of cavities fields, similar to definition 
as follows 



U) 



1 N 1 

2—1 1 



jec(i) 



(22) 



where 1 is the 2x2 identity matrix. For the physical 
fields we get 



Mi 



= tim(w \Jo, itf) + e u m (/j£ ij, u 



eeca 



iedi 



(21) 



JJ dfiw ch (v 

i=i 

J dvw ch (v)YPk J 



E 



(fc> 



A— 1 



which lead us to the self-consistency equations 



fe-i 



5 /lx-^u(i^|J ,1) -^]u(/Lt^|J,U^ 



fe>0 



8 M- u(i/|J , 1) -^u(^|J,U^) 



(23) 



Similarly, we get the following equation for the density 
of physical fields 



J 



Wphys(M) 



J]_dvew ch (ve) 



E 

fe>0 



Pk 



2 

fi -^u(i^|J , 1) -^]u(^|J,U^ 

=1 



(24) 



A. Parameter-independent equations and 
bifurcation analysis 

The above set of self-consistency eqs. (1231) do not allow 
a general analytical treatment. However, we can extract 



some information exactly. To do so, it is convenient to 
switch to the parameter-independent representation. Af- 
ter using the transformations (|18|l we obtain 
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k=0 v 1 J 

P{a) 



x 5i 



1=1 
1 



n{^ h }w ch [{p c £ h } 

2 



k-l 



H{dP* r }w gT [{pt r }}dxj e D(ij e ) 



fc-i 



k-l 



exp j /3 Jo J]] cr ■ <T£ + [3 J er • 



oo „ „ 

Wch[{P}] = 5>fc / {rfPch}W ch [{Pch}] / 
fe=0 ^ "* 



n{< r }w gr [{p g £ r }]du^(u,) 



x <5, 



(f) 



^ch JS JS 



n 



exp I /3 Jq<t • er' + /3 J cr • Uf 



and similarly for the density of physical distributions 



(25) 



OO „ 

VW{P}]=5>* / 

k=0 J 



f[{dPi}W ch [{Pi}} 



.1=1 



Y[{dP^}W gr [{P^}]d\J e D(U e 



x S 



(F) 



fldatPiiat) 



i=i 



u=i 



exp /3 J cr ■ ct£ + f3 J cr ■ n 



i=i 



6=1 



Even though eqs. i|25|l and (|26|) still look rather nontriv- 
ial they allow us to locate the continuous transitions in 
a clean way. Let us start by noticing that from the def- 
initions of the magnetization and the spin-glass matrix 
overlap, the paramagnetic solution {M a = Q a (, = 0) is 
simply given by 



W phys [{P}] - 6 {F) 
which implies 

W gr / ch [{P}] = S {F) 



P(cr)- 



1 

2^ 



P(<x)- 



2tt 



(27) 



(28) 



(26) 



r 



To find bifurcations away from (|28|l we follow |28j 
and make the so-called Guzai expansion P(er) — > 
(27T)- 1 + A(cr), with >V gr/ch [{P}] -+ W gT/ch [{A}} and 
where normalization requires that W gr / c h[{A}] = if 
J s da A(cr) ^ 0. We then examine quantities such as 

/ s {dA}Wg r/ch [{A}]A(eri) • • • A(er r ) for r = 1, 2, and as- 
sume that close to a continuous phase transition they 
are of order 0(e r ) with e arbitrarily small. Inspecting 
the possible bifurcations in first and second order in e wc 
find the following two transition lines written in terms of 
the moments (k) and (k 2 ) of the graph degree distribu- 
tion pk (for details see similar calculations in j2^): 



P-F: 



P-SG: 



1 = cos(o-i) 



(k 2 ) (k) 



(k) 



2(k) 



h(PJo) 



1 = 



7 2 (/3J) 



(k 2 ) ~ (fc) 
(k) 



2(k) 



ig(PJa) - li(j3J Q ) 



(29) 
(30) 



where we have used the following measure for D(XJ) 



with 



D(V) 



du> K(lu) S 



U 



cos u! sin u> 
— sin lu cos u> 



(31) 



K(uj) = ^S[uj -57] + -5[qj+lu} (32) 
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for some ui £ [0, 2ir]. 

B. Parameter-dependent equations, population 
dynamics and simulations 

To obtain the behavior of the order parameters with 
respect to the control parameters (such as temperature, 
average connectivity and so on) we must solve the self- 
consistency equations numerically. In order to do so, 
we switch to the parameter-dependent representation for 
self-consistency equations ()25[) and apply population dy- 
namics. 

To do the integrals in (|23[1 it is convenient to parametrize 
the spin vector cr = (cos(</>), sin(0)) with <f> S [0, 2ir). 
Then for the cavity and physical distribution we have 
that P(c) — > P(<fi) and since P(4>) is a periodic function 
of (j) we choose a representation in Fourier modes 

P((j>\fi) ~ exp I a s cos(s0) + b s sin(s^) , (33) 

\s>l S>1 J 

with (i = (oi, a2, . . . ; 6i, &2) • ■ -)j which implies that we 
have chosen, rather than polynomials <^„(<x), plane- waves 
ip n {cr) — > (sin(n(/>), cos(n0)) with weight w(cr) = 1/tt. 
After all this has been introduced into 12- il) we use pop- 
ulation dynamics to solve them. Needless to say that 
to solve the equations for an infinite number of cavity 
fields is unfeasible: we must truncate the Fourier se- 
ries in the parameter representation l|33|l . Also, we de- 
scribe the macroscopic state of the system in terms of 

just two order parameters m = M£ + and q = 

(1/2)(Q XX + Q yy ). Numerical solutions of the eqs. Ij23(l 
using population dynamics considering only two parame- 
ters /j, = (ai , 0, . . . ; bi , 0, . . .) are shown in Fig. ^ To see 
the crudeness of the truncation we have compared the 
numerical results with Monte Carlo simulations. These 
have been done for system sizes of N = 40 • 10 3 spins 
and thermalizing with the Fast Linear algorithm |33| . In 
Fig. ^the comparison between population dynamics and 
Monte Carlo simulations is fairly good (Jo = 1/2 and 
,/ = 1), supporting the truncation scheme. 

C. Equations for a given instance and proliferation 
of pure states 

So far we have been tacitly assuming that there is only 
one pure state, or equivalently one set of cavity fields fx^ 
Vi = 1, . . . ,N and Vj e di, solving the message passing 
equations (|20() . Thus iteration of these equations does 
not depend on initial conditions. The proliferation of so- 
lutions to the equations H2(J|) signals ergodicity breaking. 
To examine ergodicity brea king we follow a method sim- 
ilar to the one presented in [34( . For a given graph real- 
ization we take two different populations of cavity fields 

(i) 1 (?) 2 

and fjb~ and let them evolve according to the 




FIG. 1: Comparison between population dynamics and Monte 
Carlo simulations for the XY small- world systems with Pois- 
son degree distribution and angular disorder K(cu) with aJ = 
jiv. We plot the magnetization m = v/ Af J + My (solid lines) 
and overlap q = \ {Q XX + Q yy ) (dashed lines) against the in- 
verse connectivity c _1 . Top lines refer to T = 0.3 while lower 
ones to T = 0.5. Markers correspond to simulation results of 
system size N = 40, 000. The critical parameter values are 
in good agreement with those obtained from our bifurcation 
analysis. 

message passing equations (1201 . If there were a single 
pure state then the system would be ergodic, and so after 
some transient, the state of the system would be indepen- 
dent of initial conditions. On the other hand, dependence 
on initial conditions for long times suggests that there are 
long range correlations contradicting the assumption of 
a single pure state (i.e. the clustering property is contra- 
vened) 35] . The difference between the two populations 
of fields is measured by introducing the following param- 
eter 

1 N 

A = ^E(IK (1) H-IK (2) H)' ( 34 ) 

t=l 

(12). 

where ' is the local magnetization for the two copies 
of the system. Thus, if A ^ the solution of the iterative 
message passing equations (|20|) is non-ergodic. Notice 
that A is invariant under global rotations between the 
two copies and therefore does not take into account the 
rotational degeneracy of the Hamiltonian. 

D. Phase diagrams 

On the top panel of Fig. [3]we find the simplest possible 
scenario of our model © with p/, = e~ c c k jk\ and ZJ = 0. 
We plot the P (paramagnetic) — > F (ferromagnetic) line 
(|29H for Jq = 0, 0.5 (with Jo = returning the results of 
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T/J 2 




SG 



T/J 2 



SG 



T/J 



FIG. 2: Phase diagrams for XY small- world networks with 
bimodal distribution K(lu) = \5[u ~ 57] + kS[u) + 57] on Pois- 



son graphs for which 



= /fc! and = (top), a> = W 



(middle) and U = \n (bottom). In all diagrams, solid lines 
represent the P— >F transition whereas dashed ones represent 
P^ SG ones (lower lines: Jo = and upper: Jo = 0.5). In 
the middle panel, where for Jo = all three phases meet, dot- 
ted lines represent the F— >SG transition and dot-dashed the 
boundary above which A > 1341 1 and non-trivial ergodicity 
breaking takes place. 



FIG. 3: Phase diagrams for XY small-world networks with 
bimodal distribution K(cu) — \S[oj — 57] + \5[uj + 57] on scale- 
free graphs for which pt ~ fc -7 with k = m, . . . , oo and Zu = 



(top), 



in (middle) and u> — ^n (bottom). We have 



taken the value of the lower cutoff of pk to be m = 1. Solid 
lines represent P-^F transition whereas dashed ones represent 
P^ SG ones (lower lines: Jo = and upper: Jo = 0.5) and 
the dotted line represents F— > SG. Small world effects play 
a significant role on the phase diagram as even for small Jo 
values the ferromagnetic area increases substantially. 



[28|). In this simple case, we clearly see the small- world 
effect: even for connectivity values c<C 1 synchronization 
is possible and the F region in parameter space becomes 
significantly enlarged. The impact of rewiring is such 
that any long range connectivity value, however small 
(but finite), can lead to coherence at a finite temperature. 
A similar result based on Monte Carlo simulations has 
been obtained in [3fll3^ |. Indicatively, even for c = 0.01 
one can find a P— »F transition temperature from l|29|l at 



a finite value (around T/J w 0.02). 

For ZZJ = hir, because of the anti-aligning forces im- 
posed by the angular disorder, the P^SG (spin-glass) 
transition always occurs first (i.e. for higher temperatures 
than the P— >F) and synchronization is impossible. The 
only possible phases are paramagnetic with m = q = 
and spin glass with m = and q ^ 0. For intermediate 
values, e.g. 57 = W we see that one can have all three 
regions in the phase diagram, although the location of 
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the F— >SG line cannot be obtained from our bifurcation 
analysis. We have assumed that this is given by the line 
segment parallel to the T/J axis connecting the triple 
point where all phases meet and T/J — 38]. This 
assumption is based on physical grounds (absence of rc- 
entrance phenomena). For 57 = jir we also see that in- 
creasing the short-range bond Jo leads very quickly to an 
effective elimination of the SG area (for this example no 
SG phase appears for at least c < 0.01). Stability within 
the ordered phase, for parameter regions where all three 
phases can appear, can be found using our methodology 
of section flV CI In the middle panel of figure [3 we show 
the boundary (dotted-dashed line) separating regions of 
A = from A > 0. Above the dotted-dashed of figure 
121 we have A > and therefore initial conditions in the 
solution obtained via iteration of l|20|) are important, sig- 
naling the presence of ergodicity breaking. 

In figure [3] we plot phase diagrams for scale- free archi- 
tectures where pk ~ fc~ 7 with m < k < oo. We have 
plotted them in the (7, Tj J)-plane for different values of 
Jo, 57. As it has been noted in for the case of Ising 
spins, we find that for exponent values 7 < 3 the system 
is always in a non-paramagnetic state and this is due to 
the divergence of the second moment (fc 2 ) at 7 = 3. For 
spin systems with Jo = one finds that a lower cutoff 
m = 1 leads to a critical "f c at zero temperature which 
ceases to exist for m > 1. However, for systems on small 
worlds (Jo 7^ 0) we see that m = 1 can no longer lead to 
such a transition as even an infinitesimally small number 
of long range short-cuts is sufficient to guarantee a fer- 
romagnetic state at zero temperature. Again, the small 
world effect is striking. In fig. |3| we compare the phase 
diagrams for Jo = and Jo = 0.5. We see that even a rel- 
atively small coupling strength Jo guarantees coherence 
at finite temperatures and also for small values of 7. As 
in the case of Poisson random graphs, without angular 
disorder (57 = 0) there can only be a P^F transition, 
whereas for 57 = ^ir only a P^SG. Notice also that by 
increasing the lower cutoff value m the area where co- 
herence occurs is further enlarged, as expected, since the 
mean connectivity on the graph also grows. The dotted- 
line separating the SG and ferromagnetic phase is based 
on the Parisi- Toulouse hypothesis [38| . Due to numerical 
difficulties with population dynamics at low 7 we were 
unable to calculate the ergodicity breaking line with suf- 
ficient accuracy in this case. 

V. CONCLUSIONS 

Complex networks constitute an important area of re- 
search by virtue of the architecture through which mi- 



croscopic elements interact with one another. Typically, 
only a vanishingly small fraction of the system, contain- 
ing a finite number of other nodes, can communicate di- 
rectly with any given node. Cooperativity on such net- 
works can vary significantly depending upon parameters 
such as the degree distribution or a superimposed small- 
world structure. As such interaction structures are com- 
monly found in a wide range of real networks in recent 
years they have attracted the attention of both experi- 
mentalists and theorists alike. From a theoretical stand- 
point the notion of 'local neighborhoods' opens up new 
directions towards theories of finite dimensional systems. 

In this paper we have adapted the techniques devel- 
oped in [28j for the cavity method to study small-world 
effects in spin systems with real-valued variables defined 
on the unit circle. Such systems are not only signifi- 
cant from a physical point of view, but are also of spe- 
cial relevance in the context of e.g. Josephson junctions 
[23l IziL l2fij and coupled oscillator systems |2fj. In special 
cases where the degree distribution is given exactly (e.g. 
Poisson or scale-free) we have performed a bifurcation 
analysis and derived phase diagrams. From our results 
we see that superimposing a small-world structure sig- 
nificantly enlarges the region in parameter space where 
synchronization is possible; even for very small short- 
range coupling strengths and in the presence of disorder 
the system is able to synchronize. For scale-free archi- 
tectures similar small-world effects arc visible again for 
very small connectivity values. Interestingly, for a lower 
cutoff of the scale free distribution m = 1 the system 
always finds itself in a synchronized state along the zero 
temperature axis. This effect is solely due to the small 
world architecture. 
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